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SmunafT  This  paper  describes  methods  for  io/erriog  msthesuidcaUy  unique  local  distributioos  of  primary  cortical  current  that  uoderty 

changes  in  the  average  pattern  of  power  of  the  ongoing  (“spontaneous”)  extracranial  magnetic  fietd  of  the  brain.  In  previous  work  we 
demonstrated  that  nuthematically  unique  solutions  to  the  inverse  problem  are  possible  for  current  sources  of  the  brain's  field,  without  assuming 
a  wwaii  set  of  current  dipoles  as  a  source  model.  In  principle,  it  is  possible  to  locate  and  delineate  patterns  of  current  of  any  configuration.  In 
practice  diis  approach  applies  to  synchronized  neuronal  activiiy,  e.g.,  activity  which  is  known  to  underiy  average  evoked  or  event-related  brain 
responses.  This  paper  extends  that  approach  to  local  changes  in  incoherent  activity,  e.g^  activity  yielding  fields  or  potentiaU  that  tend  to  be 
self-canceling  wim  averaged  over  time.  This  includes  the  spontaneous  brain  activity  normally  treated  as  background  noise  when  it  accompanies 
event-related  responses.  We  demonstrate  that  local  changes  in  this  ongoing  incoherent  activity  may  also  be  uniquely  delineated  in  space  and 
time.  The  solution  b  a  covariance  matrb  characterizing  activity  acrou  an  image  surface.  Its  diagonal  elements  represent  the  spuial  pattern  of 
mean  current  power.  Evidence  b  reviewed  indicating  that  the  dbtributioa  of  the  brain’s  magnetic  field,  due  to  both  its  synchronized  and 
incoherent  neural  activity,  b  affected  by  early  sensory-perceptual  processes  and  by  higher  cognitive  pracesses.  Hence,  in  principle,  the  ability  to 
delineate  both  kinds  of  sources  in  space  and  time  makes  it  possible  to  form  more  comprehensive  dynamic  functional  images  of  the  human  brain. 
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Recent  theoretical  developments  make  it  possible  in 
principle  to  locate  and  delineate  the  pattern  of  electri* 
cal  activity  of  the  cortex  underiying  event-related 
changes  in  the  brain’s  magnetic  field.  Given  prior 
knowledge  of  the  underlying  gecMnetry  of  the  cortex, 
solutions  to  this  problem  of  source  identification  can 
be  mathematically  unique,  despite  the  widely  held  as¬ 
sumption  that  such  solutions  must  inevitably  be  am¬ 
biguous.  As  a  practical  matter,  these  mathematically 
unique  solutions  to  the  inverse  problems  are  relevant 
to  sources  of  stable  field  patterns,  i.e.,  those  that  can 
be  recovered  firmn  background  noise  by  means  of  sig¬ 
nal  averaging.  Thus,  they  apply  to  fields  of  evoked  or 
event-related  responses.  Tlie  spatial  distribution  of 
fields  due  to  the  intrinsic  activity  of  the  brain,  i.e., 
those  that  accompany  the  event-related  fields,  change 
with  time,  and  are  therefore  considered  to  be  incoher¬ 
ent.  Owing  to  this  spatio-temporal  incoherence,  they 
tend  to  be  seif-canceling  when  averaged.  This  self- 
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cancellation  permits  the  event-related  response  to 
emerge  during  averaging.  However,  recent  experimen¬ 
tal  work  demonstrated  that,  despite  its  incoherence, 
the  level  of  the  intrinsic  activity  of  the  brain  changes  as 
subjects  perform  different  mental  tasks.  Further,  these 
task-relevant  changes  arise  from  circumscribed  regions 
of  the  brain,  and  the  regions  involved  vary  depending 
upon  the  nature  of  the  task. 

Although  complicated  and  incoherent,  the  specific 
spatial  distribution  of  the  extracranial  spontaneous  field 
is  strongly  dependent  upon  the  undert^g  geometry  of 
the  cortex  as  well  as  the  relative  amounts  of  activity 
across  its  surface.  However,  when  the  variance  (power) 
of  the  spontaneous  extracranial  field  is  measured  to 
yield  a  spatial  map  of  average  field  power,  stable 
spatial  patterns  may  emerge.  These  patterns  are  known 
to  be  dependent  upon  both  the  underlying  geometry 
and  the  relative  amounts  of  activity  across  the  rorticai 
surface.  In  this  paper  we  demonstrate  for  the  first  time 
that  mathematically  unique  solutions  to  the  inverse 
problem  also  exist  for  source  configurations  that  un¬ 
deriy  the  external  distribution  of  field  power.  The 
solution  to  this  problem  also  requires  prior  knowledge 
of  the  gemnetry  of  the  cortical  surface.  However,  math¬ 
ematically  the  proof  is  quite  different  from  the  one 
establishing  that  sources  of  event-related  fields  may 
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also  be  identified.  Solutions  to  the  invene  problem  for 
event-related  field  sources  and  for  differential  levels  of 
intfinsic  brain  activity  provide  complementary  informa¬ 
tion.  Together  they  are  capable  of  providing  a  more 
complete  picture  of  brain  activity  than  can  be  obtained 
from  either  type  of  activity  alone.  We  now  provide  a 
detailed  and  documented  account  of  these  concepts. 

Ionic  currents  flowing  within  the  brain’s  neurons  are 
accompanied  by  superiinposable  magnetic  fields  that 
encircle  entire  populations  of  concurrently  active  neu¬ 
rons.  Where  the  coherent  (synchronized)  activity  of  the 
neurons  is  limited  in  area  —  so  that  the  extent  of 
cortex  involved  is  small  relative  to  the  distance  at 
which  its  field  is  measured  —  the  external  field  is 
essentially  indistinguishable  firom  one  produced  by  a 
current  dipole.  Field  measurements  make  it  possible  to 
determine  the  3-dimensional  location,  orientation  and 
strength  of  this  equivalent  current  dipole  (Cuffin  and 
Cohen  1977;  Williamson  and  Kaufman  1981,  1987). 
Typically,  a  statistical  method  is  used  to  fit  the  ob¬ 
served  field  pattern  to  one  that  would  be  produced  by 
a  hypothetical  underlying  current  dipole.  However,  a 
dipolar  source  deduced  from  an  observed  field  is  but 
one  of  many  possible  solutions  to  the  so-called  inverse 
problem  which,  in  general,  has  no  unique  solution. 
Many  different  source  coofiguratioos  could  produce 
virtually  the  same  dipolar  field  pattern. 

Generally,  the  inherent  ambiguity  of  inverse  solu¬ 
tions  makes  it  impossible  to  be  certain  of  the  actual 
generator  of  an  observed  field  pattern,  as  the  field  may 
be  due  to  the  activity  of  an  assemblage  of  dipoles 
distributed  throughout  the  intracranial  space.  In  fact, 
Scherg  and  Von  Cramon  (1985)  and  Scherg  (1990) 
devised  a  spatio-temporal  method  in  which  best  fits  to 
a  multiple  dipole  source  is  sought.  Mosher  et  al.  (1992) 
demonstrated  that  the  linear  moment  parameters  can 
be  separated  from  the  non-Iioear  location  parameters 
using  similar  approaches.  However,  it  is  not  necessary 
to  restrict  solutions  to  those  predicated  on  specific 
source  models,  e.g.,  dipoles.  For  example,  in  recent 
work  on  magnetic  source  imaging  it  was  assumed  that 
the  current  to  be  imaged  Sowed  parallel  to  a  plane  or 
to  each  of  several  parallel  planes,  but  no  particular 
configuration  was  assumed.  In  some  of  these  studies 
(Dallas  1985;  Kullmann  and  Dallas  1987)  the  shape  of 
an  arbitrary  planar  current  distribution  producing  an 
observed  external  field  was  computed  from  the  inverse 
of  the  Fourier  transform  of  the  field.  However,  this 
approach  does  not  provide  a  unique  solution  to  the 
inverse  problem,  as  magnetic  field  measurements  alone 
do  not  contain  enough  information  to  assure  unique¬ 
ness.  Nevertheless,  within  constraints  due  to  limited 
bandwidth  of  the  spatial  frequency  content  of  the  field 
(because  of  the  distance  between  the  source  and  the 
sensor),  and  ern>n  due  to  additive  noise,  the  method 
does  yield  a  cunent  distribution  in  any  arbitrarily  se¬ 


lected  plane  or  subspace  within  a  volume.  What  is 
lacking  is  unambigudus  knowledge  of  the  space  actu¬ 
ally  occupied  by  the  primary  current  distribution  within 
the  volume.  Other  authors  describe  similar  approaches 
(Dallas  ct  al.  1987;  Kullmann  and  Dallas  1987;  Singh  et 
aj.  1984;  Roth  et  al.  1989;  Tan  et  al.  1990). 

Another  related  approach  is  that  of  Himalainen 
and  Ilmoniemi  (1984,  1991),  Kullmann  ct  al.  (1989), 
Okada  and  Huang  (1990),  and  Graumann  (1991),  who 
used  linear  estimation  methods  to  determine  the  cur¬ 
rent  distributions  from  their  field  patterns.  We  should 
also  note  the  recent  work  of  George  et  al.  (1991), 
Grcenblatt  (1991),  Okada  et  al.  (1991),  Robinson 
(1991),  Szinger  and  Kuc  (1991),  etc.  who  employed 
related  methods  to  form  magnetic  source  images.  Kull¬ 
mann  (1991),  summarizing  much  of  the  work  based  on 
this  approach,  describes  how  the  linear  estimation 
method  (within  the  bounds  of  precision  due  to  noise 
and  spatial  filtering  of  the  observed  field)  reveals  rec¬ 
ognizable  images  of  the  distributions  when  they  are 
projected  onto  an  arbitrary  plane  within  the  conducting 
volume.  However,  if  the  surfaces  onto  which  the  solu¬ 
tions  are  projected  are  truly  arbitrary,  then  the  result¬ 
ing  inverse  solution  is  not  unique.  Similar  problems  are 
associated  with  probabilistic  approaches  described  by 
Garke  et  ah  (1989),  and  by  loannides  et  al.  (1989, 
1990),  although  in  principle  they  too  can  deal  with 
current  configuratioos  on  a  plane  or  stacks  of  planes 
within  a  conducting  volume. 

It  is  obvious  that  the  actual  surface  occupied  by  the 
primary  current  configuration  is  knowable.  In  fact,  the 
major  sources  contributing  to  the  magnetoencepbalo- 
gram  are  known  to  exist  in  the  cerebral  cortex.  Fur¬ 
thermore,  ample  electrophysiologicai  evidence  sup¬ 
ports  the  assumption  that  the  primary  current  flow  is, 
on  average,  normal  to  the  surface  and  not  parallel  to  it 
(George  et  al.  1991;  Kaufinan  et  al.  1991;  Wang  et  al. 
1992).  This  a  priori  information  is  one  vital  ingredient 
of  methods  for  finding  mathematically  unique  solutions 
to  the  inverse  problem. 

The  geometry  of  the  surface  containing  the  primary 
sources  needed  to  achieve  a  mathematically  unique 
solution  to  the  inverse  problem  may  be  obtained  from 
higb-resolutioo  MRl  scans.  Moreover,  basic  electro- 
physiology  teaches  us  that  postsynaptic  potentials  lead¬ 
ing  to  axial  intracellular  current  flow,  e.g.,  in  pyramidal 
cells  within  cortical  macrocolumns,  is  predominantly 
normal  to  the  surface  of  the  cortex.  These  axial 
(primary)  currents  make  the  major  coDtribution  to  the 
field  measured  normal  to  the  surface  of  the  scalp. 
Thus,  according  to  Wang  et  ai,  (1992),  any  arbitrary 
pattern  of  current  would  consist  of  elements  of  current 
flowing  in  the  same  direction  normal  to  the  surface  of 
some  region  of  cortex.  The  problem  is  thus  constrained 
to  discovering  the  net  strength  of  this  current  flow  and 
the  area  of  cortex  it  occupies. 
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The  linear  estimation  method  of  Hamalaincn  and 
Ilmoniemi  (1984.  1991)  and  Crowley  ct  al.  (1989)  em¬ 
ploys  a  so-called  minimum-norm  criterion,  which  refers 
to  the  fact  that  the  accepted  solution  is  one  in  which 
the  net  squared  current  (the  square  of  the  current 
integrated  over  the  surface)  predicting  the  observed 
field  is  a  minimum,  as  compared  with  the  net  squared 
current  of  all  possible  candidate  source  configurations. 
It  should  be  noted  that  Hamalainen  and  Ilmoniemi 
(1984)  may  have  been  the  first  to  recognize  the  impor¬ 
tance  of  prior  knowledge  of  the  underlying  source 
geometry  However,  even  when  supplemented  with 
prior  knowledge  of  the  source  surface,  the  minimum- 
norm  criterion  they  and  others  have  employed  is  still 
insufficient  to  produtx  a  solution  that  qualifies  as 
mathematically  unique.  Wang  et  al.  (1992)  recognized 
that  the  minimizatioa  of  the  sum  of  the  squares  of  the 
differences  between  the  observed  or  measured  field 
pattern  external  to  the  scalp  and  a  theoretical  pattern 
computed  from  the  source  configuratioo  meeting  the 
minimum-norm  criterion  leads  to  a  unique  solution. 
This  minimum-norm  least-squares  (MNLS)  criterion 
falls  naturally  out  of  the  specific  type  of  the  general¬ 
ized  inverse  Wang  et  al.  adopted  for  this  problem 
(Penrose  1955).  The  earlier  investigators  who  adopted 
this  same  generalized  inverse  did  not  seem  to  recog¬ 
nize  its  least-squares  feature  and  that  it  makes  a  unique 
solution  possible. 


MNLS  solution  for  coherent  sources 

Review  of  bask  concepts 

The  MNLS  approach  can  be  simply  explained  by 
considering  a  geometry  illustrated  in  Fig.  1  which  shows 
a  source  surface  composed  of  two  planar  surfaces  at 
right  angles  to  each  other,  and  an  observation  surface 
above  and  at  right  angles  to  the  source  surface.  Since 
the  source  surface  is  presumed  to  be  known,  it  coin¬ 
cides  with  the  image  surface  on  which  the  inferred 
current  pattern,  flowing  normal  to  it,  is  to  be  located. 
Fig.  2  illustrates  a  field  pattern  on  the  observation 
surface  generated  by  single  dipole  on  the  source  sur¬ 
face.  The  image  of  the  source  computed  from  the  field 
is  shown  on  the  image  surface.  Fig.  3  shows  a  region  of 
one  plane  of  the  source  surface  containing  a  crescent¬ 
shaped  pattern  of  current  elements.  These  cunents  are 
all  flowing  in  the  same  directions  at  once.  Thus,  we  say 
that  they  are  coherent  or  synchronized.  The  plots  in 
the  middle  panels  of  Fig.  3  show  the  inveise  solutions 
across  the  image  surface.  The  inverse  solution  clearly 
resembles  the  original  source.  Wang  et  al  (1992)  pro¬ 
vide  several  examples  of  this  type  of  result,  thus  provid¬ 
ing  demonstrations  of  unique  inverse  solutions  given 
the  conditions  specified  abt^e. 

Thus  far  we  have  said  little  about  the  actual  meas- 
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Fig.  1.  Model  configuration  in  wtaicfa  the  component  of  the  magnetic 
field  Bf  normal  to  the  obstnxuion  pUmt  is  measured  by  "point 
sensors"  at  the  nodes  of  a  grid  composed  of  cells  of  sue  s.  The 
center  of  the  observation  plane  lies  at  *  «  +  2  cm.  The  source  space 
is  an  I>tbaped  surface  formed  by  two  vertical  walls.  The  top  edges  of 
the  walls  art  1  cm  below  the  obscrvatioa  plane  and  the  comer 
joining  the  walls  coincides  with  the  z-ass.  The  tmngr  surface  com 
cides  with  the  source  surface  and  inverse  compuutions  arc  earned 
out  to  deduce  the  source  denstty  on  a  grid  of  size  a. 


urements  made  on  the  observation  surface.  At  every 
position  where  a  measurement  of  the  field  is  made,  the 
actual  measure  depends  upon  the  contribution  of  each 
current  element  of  the  source,  which  in  turn  is  depen¬ 
dent  upon  the  strength  of  the  element,  its  orientation 
relative  to  the  sensing  coil  and  its  distance  from  the 
coil.  Further,  it  depends  upon  the  geometry  of  the  coil 
itself  and  the  orientation  of  the  coil  with  respect  to  the 


Fig.  2.  a;  the  field  pattern  given  by  a  single  current  dipole  of  strength 
Q>1  nA-m  normal  to  the  xz  plane  of  the  source  surface  and 
located  at  the  (a,  z)»(2,  -3)  cm.  b:  the  inverse  solution  found  on 
the  image  surface.  Each  current  dipole  represents  the  deduced 
coatributkw  of  a  cell  of  area  a^.  c:  the  saom  invetse  solutmn  is 
represented  by  isocurrent  contours,  t.e.,  contours  of  constant  cur¬ 
rent-dipole  moment  density.  The  values  are  indicated  in  units  of 
isA/a. 


Page  4 


4 


obscrvatron  surface.  The  effects  of  all  of  these  facton 
may  be  taken  into  account  through  the  concept  of  coil’s 
lead  field  (see  below).  Points  on  the  image  surface  are 
taken  to  represent  elements  of  a  curr^t  configuration 
that  would  produce  the  values  of  Bf  measured  at 
positions  r],  i  - 1,...,  m,  on  the  observation  surface, 
where  the  measures  are  weighted  by  the  sensitivity  of 
the  respective  detection  coil  at  the  various  positions. 
The  distribution  of  the  current  acrc^  the  cortex  can  be 
represented  as  an  array  of  closely  spaced  current 
dipoles  of  suitable  strength  Qj  at  the  positions  rj  on 
the  source  surface,  the  values  of  Oj,  where  j  -  U...  n 
determine  the  curved  field  pattern.  The  matrix  relat¬ 
ing  the  values  B;  at  Fj  to  a  current  dipole  of  unit 
strength  (iQjl  « 1  A-m)  at  T  is  referred  to  as  the 
lead  field  matrix  L.  A  more  technical  account  of  the 
lead  field  concept  and  that  of  the  nuinimum-norm 
least-squares  (N^LS)  approach  follows.  The  non- 
mathematical  reader  may  skip  the  equations  in  the 
sections  with  a  leading  asterisk,  but  much  of  the  text  is 
understandable  without  them. 

Lead  field  analysis.  In  common  with  ail  of  the 
other  linear  estimation  and  related  methods,  the  ap¬ 
proach  of  Wang  et  al.  (1992)  assumes  that  the  primary 
current  underlying  the  observed  field  is  intracellular 
and  that  the  secondary  (volume)  cunents  make  a  negli¬ 
gible  contnbution.  Further,  the  primary  current  is  as¬ 
sumed  to  be  distributed  within  a  finite  space  fl.  The 


law  of  Biot  and  Savart  relates  the  magnetic  induction 
B(t)  to  the  current  density  itT  within  Q 


J(r)x(?-r) 


df'* 


(1) 


In  reality,  the  number  of  observations  of  magnetic  field 
is  finite.  While  the  brain  anatomy  is  complicated,  as 
pointed  out  earlier,  in  principle,  it  has  a  known  geo¬ 
metrical  configuration  h.  It  follows  that  the  observa¬ 
tions  of  the  normal  components  of  field  Bj  are  linearly 
related  to  the  current  source  density  by 


B,- /‘Li(^)•^(^)dr'^  i-1 . m.  12) 

where  T  is  the  position  vector  wthin  the  source  region 
/2,  GTe/J).  The  vector  form  LjCD  is  known  as  the 
lead  field  (Plonsey  ^972;  Williamson  and  Kaufman 
1981).  The  lead  field  L,(r)  accounts  for  the  sensitivity 
to  the  magnetic  flux  of  the  ith  pickup  coil  at  position  f, 
produced  by  unit  source  current  density  at  T.  The  lead 
Geld  within  O  is  determined  by  the  geometry  anc' 
orientation  of  the  detecnon  coil.  For  the  present  we 
make  the  simplifying  assumptioa  that  the  detection  coil 
is  so  small  that  it  may  be  considered  to  be  a  point. 

To  prove  the  principles  presented  here,  it  is  not 
necessary  to  employ  the  actual  convoluted  cerebral 
cortex  as  a  smuce  surface.  Instead,  in  the  interest  of 
clarity,  we  employ  as  a  source  space  /}  a  simple  folded 


Fig.  3.  a:  «a  etteaded  crescent-ihaped  tource  of  uniform  curreat-dipoie  dnitity  on  the  lource  nuface  with  upper  ed|e  at  z  •  -  IS  cm  and  lower 
at  -4.S  cm.  The  inverM  miutioo  on  the  image  tur&ce  ia  diown  on  the  ri|ht  aa  a  diitributwa  of  cumnt  dipoiet  and  a«  iiocontours  for 
curreat.d4)crie  mmoent  demiiy.  The  aoiutkm  has  a  aimiiar  creaoeat-diaped  outfine,  hut  ia  weakened  and  tomewhal  spread  at  treater  depths,  b; 
the  tame  source,  but  the  current  denshy  at  the  bottom  is  4  timee  that  at  the  top.  In  this  cue  a  field  tampUng  interval  cd  s  -  0.S  cm  was  used  and 

measurements  were  made  across  s  14cmxl4cmo(Me(vationpiane. 
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surface  layer  of  umforai  thickness  tv.  Therefore,  the 
volume  integral  is  reduced  to  a  surface  integral.  If  we 
further  divide  the  surface  into  a  grid  of  n  cells  with 
area  AS',  centered  at  position  fj  and  assume  that 
current  density  is  constant  within  each  small  area,  we 
can  replace  the  integral  of  Eq.  2  with  a  summation 
sign.  Eq.  2  becomes 

E  Li{i7)  r(f7)w4S;,  i-l . m.  (3) 

Note  that  J(r^)wdS]  bears  the  dimension  of  current 
dipole;  ampere-meter.  In  matrix  representation  Eq.  3 
can  be  written  as 

B-LQ.  (4) 

where  b  and  Q  are  column  vectors  and  L  is  a  mxn 
matrix. 

As  stated  earlier,  it  is  well  known  that  the  inverse 
solution  of  Eq.  4  is  not  unique.  Mathematically,  this  is 
because  there  may  exist  a  family  of  solutions  Q*  of  the 
following  homogeneous  equation  corresponding  to  Eq. 
4 

LQ»-0.  (5) 

We  know  that  the  solutions  of  !^.  5  may  have  an 
infinite  ntunber  of  elements.  If  Q  is  an  estimated 
source  of  Eq.  4  then  any  combination  of  Q  -i-  Q*  could 
be  the  solution  of  Eq.  4.  The  component  Q*  repre¬ 
sents  the  “magnetically  silent”  sources,  which  produce 

magnetic  field  on  the  observation  plane. 

The  unique  minimum-norm  least-squares  inverse  of 
the  field.  Despite  the  obstacle  posed  by  the  non- 
uniqueness  of  inverse  solutions,  we  claim  that  the 
source  can  be  estimated  by  minimizing  the  square  of 
the  residual  error  between  the  measurements  and  the¬ 
oretical  computations  from  the  estimated  sources.  In 
the  presence  of  noise,  we  minimize  the  weighted  least- 
squares  error  (or  the  statistic) 

«  2 
m  /  B-  —  0.  \ 

tninimum  of  T.  — ^ - -  j  (6) 

where  Bj  represents  values  of  the  field  computed  from 
the  estimated  source,  B|  the  measured  values  of  the 
field,  and  a,  the  rms  noise.  The  corresponding  matrix 
expression  of  Eq.  6  is 

minimum  of  lIL'O—b’U^-  (7) 

where  ||  II  is  the  Euclidean  norm  of  a  vector,  L'  the 
weighted  matrix  with  element  Lif/o„  and  b'  the 
weighted  vector  with  element  Bi/Ci  in  the  presence  of 
noise.  In  what  follows  we  will  drop  the  prime  in  L  and 
b  for  simplicity.  However,  it  is  understood  that  they 
refer  to  the  weighted  matrix  and  vector  whenever  noise 
is  present. 

In  statistics  one  generally  deals  with  cases  where  the 
number  of  measurements,  and  therefore  the  number  of 
equations,  is  greater  than  number  of  unknowns  (m  >  n). 


and  the  unknowns  arc  always  linearly  independent.  It 
is  only  in  such  conditions  that  the  least-squares  fit  will 
give  a  unique  estimate. 

In  reality  the  number  of  unknowns  may  exceed  the 
number  of  equations.  Also,  the  n  unknosvns  may  not 
all  be  linearly  independent  of  each  other  because  of 
noise  or  an  accumulation  of  round-off  errors.  In  this 
case  the  least-squares  criterion  will  not  lead  to  a  unique 
estimate.  That  is  to  say,  there  outy  exist  a  set  of 
solutions  which  all  fulfill  the  least-squares  criterion.  If 
we  further  select  one  solution  from  such  a  set  and 
require  that  the  square  of  the  image  current  ("power") 
integrated  over  the  image  surface  also  be  a  minimum, 
the  resulting  estimate  has  the  least  residual  error  and 
has  itself  the  minirnmn  power  amon^t  aU  least-squares 
solutions 

9 

mmimufflof  Minimum  of  ^Qf  (8) 

i  ' 

Equivalently,  we  seek  the  solution  that  minimizes 
the  following  Euclidean  norm  of  the  vector  Q 

/  »  ^  Y*/* 

minimum  of  HQ  A  •  Mtoimum  of  Eq?  •  (9) 

The  above  equation  is  the  so-called  minimum-norm 
which,  as  we  pointed  out  earlier,  together  with  Eq.  6 
(or  7)  represent  both  the  least-square  residual  error 
and  the  minimum  power  of  the  inverse  solution. 

There  exist  many  different  types  of  generalized  in¬ 
verses  in  mathematics.  This  particular  type  which  has 
the  feature  of  the  least  residual  error  of  the  measure¬ 
ments  and  the  minimum  power  of  the  source  is  known 
as  the  Moore-Penrose  inverse,  or  pseudo-inverse 
(Ben-Israel  and  Greville  1974;  Barnett  1990).  We  have 
adopted  the  terminology  of  minimum-norm  least- 
squares  (MNLS)  inverse  because  the  name  implies  its 
properties.  More  importantly,  the  MNLS  inverse  is 
proven  to  be  mathematically  unique. 

The  unique  MNLS  inverse  estimate  of  Eq.  4  is  given 
by 

0-L^b,  (10) 

where  L*  is  the  MNLS  inverse  matrix  of  L.  The 
structure  of  and  how  to  emnpute  it  are  discussed  in 
detail  in  Wang  et  aL  (1992).  Here  we  mention  only  that 
the  method  of  singular  value  decomposition  (SVD)  is 
used  in  (imputing  (Press  et  aL  1986;  Ciarlet  1991). 
Nevertheless,  we  list  a  few  properties  of  the  MNLS 
inverse  which  could  be  usefid  to  the  reader  in  actual 
computations  (Ben-Israel  and  Greville  1974;  Ciarlet 
1991). 

(L*')*-(L*)^  (lU) 

L*-(L^L)*L^  (lib) 

L*-L’’(I13‘)*  (He) 
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MNLS  solution  for  incoherent  sources 

Incoherent  versus  coherent  distributed  sources 

It  is  important  to  emphasize  that  the  MNLS  method 
discussed  thus  far  applies  only  to  coherent  (synchro¬ 
nous)  patterns  of  primary  current  and,  as  recounted 
above,  many  authors  have  dealt  with  various  facets  of 
the  problem.  Synchronized  current  elements,  where 
current  flows  normal  to  the  surface  of  the  cortex, 
provide  a  reasonable  source  model  for  spatially  coher¬ 
ent  event-related  fields  (Kaufman  et  al.  1991a).  Time 
averages  are  dominated  by  the  spatially  coherent  fields, 
while  those  of  the  incoherent  background  activity  tend 
to  be  self-canceling.  Evoked  fields  are  stable  in  space, 
while  the  spatial  distributions  of  the  fields  of  the 
incoherent  background  activity  change  with  time.  How¬ 
ever,  if  the  fields  due  to  ongoing  incoherent  activity  are 
squared  to  obtain  field  power,  and  the  resulting  spatial 
patterns  of  power  are  averaged,  then  the  configuration 
of  the  averaged  power  pattern  is  determined  by  the 
underlying  geometry  of  the  source  surface,  as  well  as 
by  the  statistics  of  the  elements  of  current  of  which  its 
source  is  composed  (Kaufinan  et  al.  1991a).  Also,  if,  on 
average,  a  region  of  the  source  surface  exhibits  either 
more  or  less  asynchronous  activity  than  its  surround¬ 
ings,  the  average  field  power  pattern  is  affected  and  is 
related  to  the  location  and  shape  of  the  region  in 
question.  This  suggests  that  it  is  possible  to  locate  and 
delineate  differentially  active  regions  of  cortex,  even 
when  that  activity  is  incoherent.  In  fact,  simple  meth¬ 
ods  previously  applied  to  dipole  localization  (William¬ 
son  and  Kaufinan  1981)  were  successfully  extended  to 
apply  (o  the  localization  of  the  ceoter-of-gravity  of  a 
region  of  incoherent  activity  that  was  either  stronger  or 
weaker  than  that  of  its  incoherently  active  surround¬ 
ings  (Kaufinan  et  al  1991a). 

Many  experiments  now  show  that  incoherent  activity 
revealed  by  power  patterns  of  extracranial  fields  are 
differentially  affected  by  the  performance  of  cognitive 
tasks.  Such  differential  effects  are  localized  to  different 
regions  of  the  scalp,  depending  on  the  nature  of  the 
cognitive  task.  These  e]q>erimeots  motivated  the  pre¬ 
sent  theoretical  study.  For  example,  Kaufman  et  al. 
(1991b)  demonstrated  that  power  of  the  spontaneous 
alpha  activity  originating  in  the  right  temporal  lobe  is 
suppressed  while  subjects  scan  memory  for  previously 
beard  musical  tones.  The  duration  of  this  suppression 
increases  linearly  with  the  size  of  the  memory  set,  as 
does  reaction  time  when  performance  is  relatively  free 
of  errors  (Sternberg  1966).  Further,  while  suppression 
is  detected  elsewhere,  its  duration  is  not  related  to  the 
time  required  to  scan  memory.  The  interhemispheric 
differences  among  these  effects  are  far  more  profound 
than  *'hose  encountered  in  event-related  potential  stud¬ 
ies  (Regan  1989).  The  N1(X)  component  evoked  by  the 
same  tones  and  detected  electrically  at  Cz  does  not 


covary  systematically  with  RT.  Further,  the  magnetic 
counterpart  to  NlOO  does  not  vary  consistently  with  set 
size  either,  and  changes  in  its  amplitude  with  set  size 
differ  dramatically  across  the  hemispheres.  NlOO  is  due 
to  the  activity  of  neurons  synchronized  by  the  stimulat¬ 
ing  event  and  may  w'-il  reflect  different  attcntional 
strategies,  but  it  dots  not  reflect  operations  on  short¬ 
term  memory.  Similarly,  Kaufinan  et  al.  (1991b)  re¬ 
viewed  an  extensive  P300  literature  on  this  subject  and 
found  that  neither  the  latency  nor  the  amplitude  of 
P300  (which  is  also  due  to  coherent  activity  of  neurons) 
varies  linearly  with  set  size.  Funber,  sometimes  there  is 
no  change  at  all  in  P300  with  set  size,  although  the  RTs 
of  the  same  subjects  increase  with  set  size. 

Kaufman  et  al.  (1990)  also  found  suppression  of 
alpha  frequencies  when  subjects  scanned  memory  for 
visual  forms.  However,  in  this  case  the  effect  was 
measured  over  the  visual  areas.  This  result  was  re¬ 
cently  confirmed  by  Cycowicz  et  al.  (1992).  Also,  when 
subjects  engage  in  verbal  tasks  in  response  to  visually 
presented  words,  the  suppression  related  to  performing 
the  task  is  not  present  over  the  visual  areas  (Kaufman 
et  al.  1989),  but  it  .is  over  the  left  frontotemporal  areas 
(Cycowicz  el  al.  1992). 

Pfurtscheller  et  al.  (1977,  1988a,b)  observed  dra¬ 
matic  changes  in  alpha  power  in  the  EEG  over  differ¬ 
ent  cortical  regions,  depending  upon  the  nature  of  the 
task.  These  effects  were  often  bilateral  and  difficult  to 
localize  from  the  EEG  data.  However,  the  reference- 
free  MEG  clearly  establishes  that  these  suppression 
effects  are  local,  and  not  merely  due  to  generalized 
arousal. 

All  of  this  serves  to  motivate  the  present  paper.  It  is 
clear  that  early  sensory-perceptual  and  some  cognitive 
processes  are  reflected  in  the  coherent  event-related 
responses.  It  is  equally  clear  that  local  changes  in 
incoherent  activity  reflect  other  cognitive  processes 
and  may  even  be  related  to  perceptual  processes  as 
well.  The  two  measures  complement  each  other,  and 
neither  alone  is  an  adequate  basis  for  building  a  cogni¬ 
tive  neuroscience.  The  MNLS  method  has  now  made  it 
possible  to  delineate  the  cortical  areas  involved  in  the 
coherent  aspect  of  brain  activity.  To  complete  the 
picture  we  now  extend  the  MNLS  method  to  make  it 
possible  to  find  unique  inverse  solutions  to  the  prob¬ 
lem  of  delineating  the  cortical  areas  exhibiting  changes 
in  incoherent  activity. 

The  mmimwn~norm  least-squares  inverse  of  power. 
The  basic  problem  attendant  upon  identifying  differen¬ 
tial  incoherently  active  regions  is  that  we  deal  with 
average  field  power  as  the  basic  measure  and  not  field 
per  se.  We  introduce  a  novel  mathematical  develop¬ 
ment  for  solving  the  inverse  problem  for  such  meas¬ 
ures. 

Previously  the  values  of  the  field  Bj  measured  nor¬ 
mal  to  the  observation  surface  at  m  positions  were 
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used  to  fonn  a  column  vector  b  of  m  elements.  If  we 
now  introduce  an  auto-correlation  matrix  Bp  in  terms 
of  the  vector  b,  we  obtain 

with  elements  (Bp)^  «  Bj  •  Bj,  where  i,  j  -  1,  2,.-,  m. 
Therefore,  the  diagonal  elements  of  Bp  are  the  field 
power  Bf,  rather  than  field. 

Manipulations  of  Eq.  4  reveal  a  new  linear  model 
which  relates  field  power  to  the  source  power,  as 
opposed  to  a  linear  relation  between  Geld  per  se  and 
source  strength  as  given  by  Eq.  4.  Thus, 

Bp-LQpL^  (13) 

where  Qp  is  the  auto-correlation  matrix  of  source 
image,  (Qp)u  -  •  Q,,  1  -  1, 2,...,  n,  with  the  diago¬ 

nals  as  the  power  of  the  source.  Note  the  difference 
between  Eq.  4,  which  has  a  matrix  and  two  vectors,  and 
Eq.  13,  which  contains  all  matrices.  Fortunately,  for 
the  matrix  equation  Eq.  13  there  also  exists  a  unique 
MNLS  inverse.  With  the  help  of  the  MNLS  inverse 
theory  (Penrose  19S5;  Ben-Israel  and  GreviUe  1974) 
and  the  property  given  in  Eq.  11a,  we  are  able  to 
derive  the  MNLS  inverse  of  the  source  power  for  a 
given  field  power  distribution 

Qp-L*Bp(l-/  (14) 

Applying  this  equation  makes  it  possible  to  delin¬ 
eate  regions  of  cortex  whose  levels  of  incoherent  activ¬ 
ity  deviate  firom  a  baseline  because  of  some  ongoing 
cognitive  process,  or  even  because  of  effects  of  tran¬ 
sient  ischemia,  cortical  hyperactivity,  abnormally  low 
metabolic  level,  etc. 

If  we  use  <>  to  denote  the  time  average,  and  take 
the  time  average  of  the  field  power  given  by  Eq.  14,  we 
obtain 

<Qp>-L*<Bp>(L*)^.  (15) 

Note  that,  the  lead  field  matrix  L  is  time  invariant, 
and  so  is  its  MNLS  matrix  L'^,  provided  that  the 
measures  made  at  different  times  are  made  at  the 

same  positions.  This  is  the  basis  for  the  extended 

MNLS  inverse. 

It  is  important  to  note  that  the  concept  of  direction 
of  current  flow  is  totally  meaningless  when  dealing  with 
field  power  related  to  such  phenomena.  These  differ 
from  evoked  responses  and  similar  time-evolud  events 
where  direction  of  underlying  current  flow  is  of  consid¬ 
erable  importance. 


Methods  and  results 
/mage  representation 

An  inverse  solution  for  a  set  of  field  measurements 
can  be  represented  in  two  convenient  ways.  One  ts  by 


an  array  of  cunent  dipoles  oriented  perpendicular  to 
the  image  surface,  each  dipole  placed  at  the  center  of  a 
cell  of  area  a^.  The  image  is  thus  represented  by  the 
distribution  of  strengths  of  the  dipoles,  each  represent¬ 
ing  the  current  dipole  moment  per  cell  (as  shown  in 
Fig.  2b).  Another  representation  is  the  pattern  of  iso¬ 
contours  on  the  image  surface  that  describe  the  dipole 
moment  density  (as  shown  in  Fig.  2c).  Both  of  these  are 
sometimes  referred  to  as  the  "current  image"  obtained 
fiom  the  data. 

An  inverse  solution  for  a  set  of  time-averaged  field 
power  measurements  can  be  represented  in  analogous 
ways.  Instead  of  an  array  of  current  dipoles  across  the 
image  surface  whose  moments  are  specified,  we  use  an 
array  whose  time-averaged  square  moments  are  speci¬ 
fied  (as  shown  in  Fig.  6b).  Arrows  directed  toward  the 
viewer  indicate  positive  strength  for  the  inverse  solu¬ 
tion.  Also,  the  image  can  be  specified  by  the  pattern  of 
isocontoun  that  describe  the  mean  square  moment 
density.  Both  of  these  can  be  referred  to  as  the  “cur¬ 
rent  power  image." 

Simulations 

The  results  of  this  theoretical  work  are  in  the  form 
of  simulations  in  which  surfaces  representing  the  cere¬ 
bral  cortex  (or  any  other  thin  layer  of  spontaneously 
active  neural  tissue  within  the  brain)  are  populated  by 
a  large  number  of  perpendicular  primary  current  ele¬ 
ments  (dipoles)  of  random  orientatioo  and  magnitudes. 
Subsets  of  these  elements  are  either  mcremented  or 
decremented  in  magnitude,  and  the  net  fields  of  all  of 
the  elements  are  summed  at  the  observation  surface 
(Fig.  1).  It  is  to  be  noted  that  the  simulations  described 
here  are  merely  illustrations  of  bow  the  extended 
MNLS  method  may  be  employed,  and  are  not  of  them¬ 
selves  proof  of  the  uniqueness  of  the  solution.  The 
uniqueness  is  inherent  in  the  mathematics  of  the  MNLS 
inverse. 

In  the  simulation,  the  current  dipole  moments  popu¬ 
lating  the  source  surface  are  sampled  from  a  uniform 
distributmn,  and  a  new  random  seed  is  applied  to  the 
random  ntunber  generator  prior  to  the  selection  of  any 
array.  This  assures  an  ever-changing  field  pattern  at 
the  observation  surface.  The  field  at  each  detector  is 
squared  to  obtain  field  power  (Eq.  12).  One  hundred 
such  plots  representing  different  time  series  are  aver¬ 
aged  together  to  form  the  average  power  plot.  Then 
the  extended  MNLS  method  is  applied  to  deduce  the 
image  of  the  source.  This  image  is  compared  with  the 
configuratioo  of  the  original  source  to  determine 
whether  the  inverse  solution  is  a  reasonable  estimate 
of  the  actual  source  which,  in  this  case,  is  either  the 
incremented  or  decremented  region  of  incoherent  ac¬ 
tivity. 

In  the  various  sunulatioos  carried  out  here,  the 
background  activity  is  spread  over  the  entire  source 
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surface  (the  L-shaped  wall  of  Fig.  1).  Each  wall  of  the 
surface  is  assumed  to  be  5  cm  x  5  cm.  The  randomized 
magnitudes  and  directions  of  cunent  flow  (to  simulate 
background  “noise")  are  given  values  that  range  from 
- 1  -*  + 1.  It  is  important  to  note  that  there  is  one 
major  difference  between  the  walls  of  our  source  model 
and  the  actual  cortex.  The  contours  of  the  latter  are 
rounded  and  not  sharply  terminated.  Rounded  con¬ 
tours  of  the  actual  cortex  insure  a  gradual  rotation  of 
its  macrocolumns,  with  those  in  a  sulcus  being  largely 
tangential  in  orientation  with  respect  to  the  surface, 
and  those  approaching  a  gyrus  tending  to  become  more 
nearly  radial  in  orientation.  Since  the  latter  produce 
much  weaker  fields  at  the  surface  than  do  the  tangen¬ 
tial  current  elements  within  the  wall  of  a  sulcus,  there 
is  a  graduated  fall-off  of  the  contnbutioos  of  elements 
nearer  the  skull.  To  simulate  this  in  our  model  we 
allow  a  smooth  attenuation  of  the  magnitudes  of  the 
stiU  randomly  selected  dipoles  beginning  about  0  J  cm 
away  from  the  top  and  side  edges  of  the  L-shaped  wall 
and  reaching  zeros  at  the  edge. 

As  indicated  above,  using  different  initial  random 
seeds,  100  samples  of  source  configuratioos  were  gen¬ 
erated.  The  bottom  of  Fig.  4  is  one  of  the  samples. 
Field  values  are  computed  for  each  sensor  on  the  12 
cm  X  12  cm  observation  plane,  which  is  1  cm  above  the 
top  edge  of  the  source  surface.  In  this  simulation  the 
sensors  are  spaced  about  1  cm  apart  As  indicated,  the 
100  plots  of  field  power  thus  generated  are  averaged. 
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Fig  4.  a:  avenge  held  power  acron  the  observatioa  plane  Crom  a 
random  array  of  dipoles  with  nrengtha  ranging  from  -1  to  -»-t 
computed  for  100  such  independent  random  samplea,  one  of  which  it 
illustrated  in  (b).  c:  contour  plot  of  the  avenge  field  power,  whoee 
two  extensions  lie  over  the  walJt  of  the  source  tptce. 


Fig.  5.  Simulation  for  avenge  field  power  computed  in  the  sane 
manner  as  for  that  of  Fig.  4,  except  that  the  underlying  source 
surface  wai  of  a  more  complicated  shape.  The  croas-Iike  shape 
shewn  below  is  the  outline  of  fissures  4  mm  wide  sepanting  "corti- 
cai”  walls  about  2  cm  long  and  2  cm  wide  and  enending  about  3  cm 
into  a  model  "bead’'  (Kauhnan  et  al.  1991a). 


This  is  Ulustiated  by  the  3-D  power  plot  in  Fig.  4a.  It  is 
of  some  interest  to  compare  this  plot  with  one  of  those 
depicted  in  Kaufman  et  al.  (1991a)  shown  here  as  Fig. 
5,  where  the  underlying  cortex  was  simulated  by  a 
cross-shaped  arrangement  of  sulci.  There  was  clearly  a 
4'Iobed  power  pattern  which  was  topologically  related 
to  the  shape  of  the  underlying  structure.  In  this  case 
the  average  power  pattern  resembles  the  shape  of  the 
underlying  l-shaped  structure  —  once  again  demon¬ 
strating  the  dependency  of  average  field  power  on 
structure.  It  is  worth  noting  in  passing  that  this  con¬ 
firms  our  earlier  conclusion  that  the  actual  geometry  of 
the  cortex  must  be  explicitly  taken  into  account  if  one 
is  to  understand  scalp-detected  phenomena,  whether 
they  be  potentials  or  fields.  In  any  event,  in  this  case 
the  average  field  power  plot  reveals  the  existence  of 
the  underlying  L-shaped  source  surface  when  the  inco¬ 
herent  background  activity  is  uniform  on  average  over 
time.  As  we  shall  see  bekrw,  departures  from  uniform¬ 
ity  also  affect  the  shape  of  the  pattern. 

Incremented  incoherent  activity 

As  a  first  test  of  the  method,  on  average,  the  current 
dipole  moment  (amplitude)  of  a  small  circular  region 
of  radius  0.6  cm  on  one  wall  is  enhanced  by  a  factor  of 
10  as  compared  with  the  otherwise  uniformly  dis¬ 
tributed  activity  of  its  surroundings,  which  is  as  de¬ 
scribed  above.  Fig.  6a  shows  one  sample  of  100  such 
source  distributions.  Note  that  the  arrows  signify  cur¬ 
rent  elements  where  the  directions  of  current  flow  are 
randomly  related  to  each  other,  unlike  the  sources  of 
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Fig.  6.  a:  the  source  surface  of  Fig.  1  is  populated  with  a  large 
number  of  current  dipoles  of  randomly  seledted  orientatioas  and 
strengths,  but  those  within  a  smaU  (0.6  cm  radhis)  regno  ate,  on 
average,  10  times  stronger  than  those  in  the  surrounding.  The 
current  power  image  distribution  (b)  and  isopower  density  contours 
(c)  are  shown  on  the  image  surface.  The  arrow  length  in  (b)  is 
proportional  to  the  average  power  at  each  location. 


evoked  responses.  We  first  computed  the  average  field 
power  distribution  at  the  sensors  on  the  observation 
surface  and  the  correlations  among  different  sensors 
(Eq.  12).  Then,  in  a  one-step  process,  we  cmnputed  the 
extended  MNLS  inverse  over  the  entire  image  surface. 
According  to  Eq.  IS,  if  the  time  average  of  field  power 
is  used,  rather  than  field,  the  resulting  inverse  solution 
gives  the  time-averaged  current  power  image. 

Since  the  extended  MNLS  approach  guarantees  that 
the  result  is  the  best  (in  the  sense  of  the  least-square 
residual  error  between  the  measurements  and  the  esti¬ 
mated  field  values)  an  iterative  search  process  in  which 
residual  errors  are  compared  at  each  step  to  find  a 
minimum  is  unnecessary  (as  in  a  typical  non-linear 
least-squares  method).  Fig.  tib  is  a  perspective  view  of 
the  current  power  image  derived  from  the  extended 
MNLS  inverse.  The  power  at  each  location  is  indicated 
by  the  length  of  an  arrow,  whose  direction  is  arbitrarily 
chosen  to  be  toward  the  reader.  Heads  are  included  to 
better  illustrate  their  length.  Fig.  6c  is  an  isocontour 
plot  representing  the  same  inverse.  Note  that  the  con¬ 


tours  are  centered  on  a  point  at  the  coordinates  (x, 
z)-(2,  -3)  cm,  which  was  the  exact  center  of  the 
original  incremented  current  distribution  of  the  source 
surface. 

We  next  consider  the  case  of  two  incremented  circu¬ 
lar  regions,  located  on  the  source  surface  at  the  coordi¬ 
nates  (x,  z)w(2,  -3)  cm,  and  (y,  z)-(2,  -3)  cm. 
While  the  current  elements  vary  at  random  over  time 
inside  as  well  as  outside  the  incremented  regions,  the 
ratio  of  averaged  amplitude  of  the  two  incremented 
regioiu  relative  to  their  sunoundings  is  10:1.  One  of 
the  100  source  samples  is  shown  in  Fig.  7b.  Fig.  7a  is 
the  plot  of  the  field  power  averaged  over  ail  100 
samples.  It  is  of  some  interest  to  compare  this  plot  with 
the  L-sbaped  plot  of  Fig.  4,  which  is  based  on  an 
average  of  lOO  samples  of  randomly  selected  current 
elements  where  there  is  no  net  increment  or  decre¬ 
ment  of  any  region.  When  two  regions  are,  on  average, 
more  “active”  than  their  surrounding,  the  L-sbaped 
plot  is  transformed  into  one  with  three  lobes,  thus 
illustrating  bow  both  the  geometry  of  the  source  sur- 
fiice  and  the  statistics  of  its  activity  affect  the  distribu- 


Fif.  7.  Simitw  to  Fig.  6,  except  that  two  0.6  cm  radius  regions,  one  on 
each  wall  of  the  louice  surface,  have  cumnt  dipole  momenu  incre- 
meoted  by  an  average  fiKtor  of  10.  The  avenge  power  distributios  is 
computed  for  100  random  arraya  of  dipole  momeoa.  The  bonom  of 
the  figure  shows  the  current  power  imnge  distributioa  determined 
from  the  tverage  field  power  pattern  and  the  known  stupe  of  the 
source  surface. 
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Fig.  8.  a:  2  of  the  100  samplea  used  io  computing  avenge  field  power  for  a  aeacent-shaped  area  of  incoberent  activity,  (b)  The  inverse  soiuiton 
representing  the  spatial  distribution  of  current  power  and  (c)  corresponding  isodensity  current  power  plot. 


tioo  of  field  power  at  the  observatioo  plane.  Fig.  7c  is 
the  inverse  solution  computed  using  the  extended 
MNLS. 

To  iiVostrate  the  potential  power  of  this  method  in 
delineating  the  shape  of  a  region  of  differential  inco* 
herent  activity,  a  crescent-shaped  portion  of  the  source 
surface  similar  to  that  of  Fig.  3  contains  current  ele¬ 
ments  which,  on  average,  are  10  times  stronger  than 
those  of  the  surrounding  region  of  the  surface.  The  top 
of  the  enhanced  region  is  at  z  -  -  2.5  cm  and  the 
bottom  at  z-  -4.5  cm.  Fig.  8a  shows  2  of  the  100 
samples  of  the  random  distributions  used  in  computing 
the  average  field  power  at  the  observation  plane.  As 
before,  each  of  these  samples  was  created  by  applying 
different  initial  seeds  to  the  random  number  generator. 
However,  the  current  power  image  (Fig.  8b)  provides  a 
reasonable  rendition  of  the  crescent  shape.  This  is 
made  clearer  in  Fig.  8c  by  the  isocontours  of  the 
current  power  image. 

While  the  enhanced  region  in  the  preceding  exam¬ 
ple  was,  on  average,  always  10  times  that  of  its  sur¬ 
roundings,  we  also  experimented  with  different  ratios. 
Thus,  the  3  plots  in  Fig.  9  are  current  power  images  in 
which  a  circular  region  0.6  cm  in  radius  is  incremented, 
but  the  amplitude  enhancement  is  5  times  (Fig.  9a),  4 
times  (9b),  and  3  times  (9c)  that  of  the  surroundings. 
Note,  these  are  all  drawn  with  the  same  scale,  while 
the  10:1  plot  of  Fig.  6  is  shown  with  a  different  scale. 

If  the  circular  region  is  moved  (k>wnward  away  from 
the  observation  plane  by  15i  cm,  the  center  of  the 
region  is  now  at  z  •  -4.5  cm  and  the  magnetic  field  at 
the  observation  plane  is  greatly  weakened.  In  this  case. 


Fig.  9.  ImMc  dmu  power  diitnbutioDi  when  a  regioa  0.6  era  in 
radhu  cemeicd  tt  t  depth  z  •  3  an  a  iaooheteady  active  with  dipole 
axmeats  of  current  elemeoB  increoented  by  a  (Ktor  of  (a)  S.  (b)  4, 
and  (c)  3  times  the  level  of  the  nurouadiiiii.  The  power  levels 
expreiied  in  the  imagea  grow  profre*ah>ety  weaker  with  the  reduced 

ratio. 
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Fig.  10.  Similar  to  Fig.  9,  eicept  that  the  sources  are  at  a  depth 
z  ••  4  J  cm  aod  dipole  moments  are  mcremeated  by  a  factor  of  (a)  7, 
(b)  6  and  (c)  5  times  the  level  of  the  surroundings. 


if  the  incretnented  region’s  strength  is  less  than  5  times 
that  of  its  surroundings,  the  incremented  region  is  no 
longer  discriminable  in  the  current  power  image  (Fig. 
10).  The  gradual  deterioration  in  the  ability  to  visualize 
current  power  images  for  incremented  regions  of  a 
particular  size  is  illustrated  in  Fig.  10  where  the  ratios, 
from  top  to  bottom,  arc  7:1,  6:1, 5:1,  respectively.  Note 
that  Figs.  9  and  10  were  plotted  at  the  same  scale. 
These  examples  are  meant  merely  to  illustrate  the  fact 
that  inverse  solutions  nmy  not  be  successful  for  regions 
with  levels  of  differential  activity  that  are  too  small 
considering  the  distance  of  the  detection  coil  from  the 
source  area,  or  when  the  area  is  smaller  than  some 
minimum  size.  The  degree  to  which  these  interacting 
parameters  will  affect  the  application  of  the  extended 
MKLS  method  in  real  situations  is  clearly  an  empirical 
question. 

Because  of  these  potential  problems  we  devised  a 
variation  on  the  basic  method  described  thus  far.  This 
variation  significantly  enhances  the  ability  to  visualize 
“target”  areas  in  current  power  images  when  signal- 
to-notse  ratio,  area,  and  distance  were  not  favorable,  as 
in  the  foregoing  examples.  In  this  variation  we  first 
evaluate  the  field  power  associated  with  the  overaU 
background  when  the  target  area  is  absent  (as  shown  in 


Fig.  4).  An  example  of  this  in  a  real  situation  would  be 
to  plot  the  power  of  the  component  of  the  field  normal 
to  the  surface  of  the  scalp  of  a  resting  human  subjeci. 
We  then  find  the  extended  MNLS  solution  for  the 
current  power  image  to  define  a  baseline  image  of  the 
background.  Then  the  localized  area  is  imposed  on  this 
background  as,  for  example,  when  the  subject  becomes 
engaged  in  a  mental  task  of  some  kind,  thus  altering 
the  level  of  activity  of  one  or  more  regions  of  the 
cortex.  We  then  find  the  extended  MNLS  for  the 
current  power  image  with  this  target  area  or  areas 
present.  Subtracting  the  image  of  the  baseline  distribu¬ 
tion  from  the  image  containing  the  target  area  gives  a 
current  power  image  difference  plot,  such  as  those 
illustrated  in  Fig.  11.  The  directions  of  the  arrows  in 
the  image  surface  reflect  the  sign  of  the  difference 
between  the  suppressed  and  baseline  activity.  Positive 
values  are  indicated  by  arrows  directed  toward  the 
reader.  Note  that  the  target  in  Fig.  11a  is  one  that  is 
difficult  to  visualize  in  Fig.  lOc.  Note  also  that  the 
noise  in  the  surroundings  is  much  improved,  even 
though  accentuated  by  the  greater  magnification  of  a 
factor  of  2.5,  the  scale  used  to  plot  the  image  strength 


Fit-  11.  laage  power-difference  ploa  for  the  sources  of  Fig.  10, 
where  the  imafe  of  the  backirouod  is  subtracted  from  the  image  of 
oackground  plus  setive  dreuiar  regioa.  for  diptrie  moaienis  in  the 
active  region  that  are  incremented  by  a  factor  of  (s>  S,  (b)  4  and  (c)  3 
times  the  level  of  the  surroundinp. 
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in  Fig.  11-  It  ‘S  important  to  emphasize  that 
different  initial  seeds  were  been  used  for  simulating 
the  baseline  activity  and  for  simulating  the  activity 
containing  the  targets.  This  indicates  that  a  similar 
approach  may  be  used  in  real  situations  where  neu¬ 
ronal  activity  is  simply  not  going  to  have  a  constant 
distribution  over  long  periods  of  time,  despite  the 
long-term  stability  of  average  MEG  and  EEG  power. 

Decremented  incoherent  activity 

Thus  far  we  have  dealt  only  with  target  areas  whose 
activity  has  a  somewhat  higher  level  than  that  of  the 
surroundinp.  However,  this  is  a  distinction  based  on 
an  assumed  convention.  It  is  equally  logical  to  consider 
these  patterns  as  those  in  which  the  so<alled  sur- 
roundinp  is  decremented  as  compared  to  the  smaller 
target  areas.  So,  if  very  large  areas  of  cortex  should 
exhibit  alpha  suppression,  for  example,  then  the  meth¬ 
ods  we  have  used  would  be  able  to  delineate  those 
areas.  Moreover,  we  may  also  consider  the  reverse 
situation,  in  which  the  suppressed  target  region  is  small 
relative  to  the  surroundings,  since  this  may  have  impor¬ 
tant  applications.  For  example,  relatively  small  cortical 
regions  affected  by  ischemia  may  exhibit  less  overall 


Fig.  12.  a:  ninpic  of  a  source  diatributioa  with  a  circular  regioa  that 
it  decremeated  b.  c:  iraage  power  dittrSnitioiia  ''btaloed  by  aubtract- 
ing  the  average  of  100  such  aaoipiea  from  the  baaeiine  inverte  of  100 
different  samplea  of  uiiifdxnn  activity,  begiiiiung  with  (hffereot  leedi. 


Rf.  13.  a;  creacent.ahape<i  aouroe  region  of  tuppreeaed  activiry.  b; 
image  power  diatrib«k»  difference  obtained  by  tubtracting  the 
current  power  image  of  100  taavica  of  (a)  from  the  current  power 
image  of  100  aampiea  of  unilbnn  activity  obtained  (ram  differeoi 
initial  seeds,  c;  ioage  power  diatrtnition  dihereooe  obtained  as  ui  (b) 
but  narnng  horn  the  same  initiai  seeds. 

activity  than  their  surroundinp.  Also,  as  implied  in  the 
work  of  Kauffnan  et  ai.  (1991b),  relatively  circum- 
senbed  areas  may  exlubit  suppression  as  subjects  en¬ 
gage  ir  scanning  short-term  memory  for  tonal  stimuli. 
It  would  be  of  considerable  value  to  locate  and  delin¬ 
eate  such  areas. 

One  of  the  difficulties  associated  with  deteaing 
small  decremented  target  areas  is  that  the  extracranial 
field  is  dominated  by  the  background  activity  which 
originates  in  a  relatively  much  larger  area  of  cortex. 
This  could  result  in  an  extremely  disadvantageous  stg- 
nal-to-noise  ratio  which  may  make  it  difficult  to  detect 
any  change  simply  by  examining  at  the  inverse  plots  of 
current  power  image.  While  in  actual  suppression  data 
studied  by  Kaufinan  et  al.  (1990, 1991a)  the  empirically 
observed  suppression  is  very  pronounced  and  may  of¬ 
ten  be  seen  in  single  trials,  this  problem  could  arise  in 
the  case  of  subtler  effects.  To  deal  with  this  in  our 
simulation,  we  begin  with  a  circular  target  area  having 
a  radius  of  1  cm.  On  average,  the  strength  of  activity 
within  this  target  is  1:10  that  of  its  surroundinp.  A 
single  sample  of  a  source  distribution  is  shown  in  Fig. 
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12a  In  our  computations,  100  such  sajnples  with  differ¬ 
ent  initial  seeds  were  used  in  the  simulation.  The 
current  power  image  when  the  source  was  decre¬ 
mented  was  subtracted  from  the  current  power  image 
of  the  baseline.  Fig.  12b  and  c  show  these  power  image 
difference  plots.  Note  again  that  different  initial  seeds 
were  used  for  producing  the  inverse  power  with  a 
decremented  region  and  the  baseline  inverse.  This 
produces  patches  of  activity  in  the  image  that  differ 
from  one  computation  to  another. 

The  above  procedures  were  employed  to  detect  a 
crescent  shaped  region  of  suppressed  activity  (Fig.  13a) 
with  the  ejctended  MNLS  inverse.  Under  the  condi¬ 
tions  of  this  simulation,  where  the  target  pattern  ex¬ 
tends  over  a  wide  range  in  depth  (with  the  top  at 
z  »  -  2.5  cm  and  the  bottom  at  z  -  -  4.5  cm)  the  lower 
part  of  the  current  power  image  is  buried  in  the  noise 
(Fig.  13b).  Of  course,  where  the  surrounding  activity  is 
more  stable  over  tune,  it  is  possible  to  obtain  an  image 
of  higher  quality.  For  example,  the  plot  in  Fig.  13c  is 
one  where  the  baseline  of  uniform  activity  and  the 
distributions  containing  the  target  shapes  were  created 
by  applying  the  same  initial  seed  to  the  random  num¬ 
ber  generator.  The  resulting  power  image  difference 
plot  shows  a  well-defined  crescent,  because  there  is  an 
“exact”  cancellation  of  the  background  activity.  The 
poorer  image  of  Fig.  13b  results  from  no  correlation  of 
baseline  and  target,  while  the  higher  quality  image 
results  from  perfect  correlation.  Empirically,  the  result 
will  probably  be  intermecuate  between  these  two  solu¬ 
tions.  In  particular,  the  length  of  time  over  which 
recording  are  made  will  have  a  substantial  influence  on 
the  quality  of  the  subtraction  technique. 


CoQclusioas 

The  main  conclusion  vf  this  paper  is  that  it  is 
possible  in  principle  to  estimate  mathematically  unique 
inverse  solutions  to  recover  the  location,  shape  and 
magnitude  of  a  dififerentially  active  region  of  cortex, 
even  when  the  activity  is  incoherent.  This,  together 
with  earlier  work,  makes  it  possible  to  uniquely  define 
both  coherent  and  incoheient  activity  using  maps  of 
the  extracranial  field  and  its  power  together  with 
MRI-based  reconstructions  of  the  individual  subject’s 
or  patient’s  brain. 

Many  problems  remain  to  be  overcome  before  the 
ide^  developed  here  can  be  employed  in  practice.  For 
example,  the  effects  of  measurement  erron  and  in 
accuracy  of  information  about  uodertying  cortical  ge¬ 
ometry  have  yet  to  be  fully  ascertained.  Our  own 
ongoing  work  suggests  that  an  accuracy  of  better  than  4 
nun  in  knowledge  of  cortical  geometry  may  be  required 
(Wang  1992;  Wang  et  al.  1992).  This  is  possible  with 
the  relatively  rer«ot  advent  of  low-distortion  1  mm 


resolution  MRI  systems.  However,  as  a  practical  mat¬ 
ter  it  is  not  yet  possible  to  automatically  segment 
cortex  from  MRI  scans  and  use  these  regions  of  inter¬ 
est  in  constructing  accurate  3-dimensional  representa¬ 
tions  of  the  cortex.  Manual  methods  are  widely  used, 
and  these  are  inherently  innaccurate.  Work  towards 
resolving  these  problems  is  underway  in  our  and  many 
other  laboratories.  Further,  effects  of  noise  must  also 
be  considered  in  more  detail.  Our  work  (Wang  1992; 
Wang  ct  al.  1992)  shows  that  the  effects  of  increased 
noise  may  be  offset  by  increasing  the  numbers  of 
positions  at  which  measurements  arc  made,  but  this 
too  must  be  investigated  further.  Extracortical  sources 
of  noise  require  special  consideration.  At  the  present 
time,  instrument  noise  is  largely  negligible,  as  modem 
ihin-fihn  SQUID  sensors  provide  noise  levels  as  good 
as  5  fTv^Hz ,  which  is  a  factor  of  4  better  than  older 
more  widely  used  neuromagnetometer  systems.  In 
dealing  with  intrinsic  activity,  the  fields  of  interest  may 
range  upwards  of  1  pT,  which  is  much  larger  than  this 
instrument  noise  level.  Biological  noise  sources,  e  g., 
heart  and  eye  movements,  are  entirely  negligible  with 
second-order  gradiometers  with  a  4-5  cm  baseline 
(Romani  et  al.  1982),  but  can  be  consequential  when 
using  long  baseline  first-order  gradiometers.  Heart  sig¬ 
nals  are  unlikely  to  be  sensed  by  planar  gradiometers, 
for  example  HamiUainen  1989).  Therefore,  interfer¬ 
ence  firom  extracortical  sources  is  a  technical  problem 
that  can  be  dealt  with  by  using  an  appropriate  pickup 
coil  design  or  by  use  of  adaptive  filtering  techniques. 
Changes  in  the  position  of  cortex  time-locked  to  pul¬ 
satile  changes  in  blood  volume  in  the  brain  are  not 
likely  to  be  of  such  a  magnitude  as  to  be  of  importance. 

Of  greater  importance  is  the  fact  that  the  adjacent 
walls  of  the  brain’s  sulci  are  often  in  close  proximity  to 
each  other,  so  that  it  may  not  be  possible  to  accurately 
determine  the  wall  to  wluch  some  change  in  level  of 
ongoing  activity  actually  belongs.  Simulation  studies  in 
which  actual  cortical  contours  are  used  may  clarify  the 
degree  to  which  this  is  a  problem,  and  also  the  degree 
to  which  various  image  enhancement  techniques  may 
be  employed  to  sharpen  the  resolution  and  bow  changes 
in  the  density  of  sampling  the  extracranial  magnetic 
field  can  affea  it.  Analogues  to  the  strategies  em¬ 
ployed  in  image  processing  are  worthy  of  further  explo¬ 
ration. 

Assuming  that  the  differential  levels  of  activity  of 
the  different  regions  of  the  brain  are  significantly  re¬ 
lated  to  sensory-perceptual  and  cognitive  processes, 
and  all  the  evidence  accumulated  thus  far  suggests  that 
they  are,  then  approach  described  in  this  paper  makes 
it  possible  to  envision  a  new  kind  of  brain  mapping. 
We  propose  to  dub  this  mapping  dynamic  functional 
imaging  because  it  will  rev^  changes  in  activity  in 
different  regions  with  a  temporal  resolution  measured 
in  fractions  of  a  second.  LK  put  milliseconds  here,  but 


Page 


14 


alpha  does  oot  change  much  in  1  ms.  PET  blood  flow 
methods  also  provide  fuo<^oaal  images,  but  the  time 
resolution  is  no  better  than  about  40  s,  and  data  must 
be  averaged  over  subjects.  MRI  techniques  show 
promise  of  better  temporal  resolution,  but  this  possibil¬ 
ity  remains  to  be  demonstrated. 

While  we  have  restricted  our  discussion  to  the  brain's 
magnetic  field,  it  may  well  become  possible  to  develop 
similar  methods  in  which  electroencephalographic  data 
are  used  instead.  In  principle  there  is  no  reason  why 
this  cannot  be  accomplished,  provided  that  an  accurate 
picture  of  the  electrical  properties,  e.g.,  conductivities, 
of  the  individual  patient’s  skull,  brain,  and  other  tis¬ 
sues  that  affect  the  fiow  of  volume  currents  are  taken 
into  account  when  computing  the  solution  for  the  in¬ 
verse  problem. 

Finally,  although  we  have  emphasized  the  human 
cerebral  cortex  throughout,  it  should  be  borne  in  mind 
that  there  is  no  reason  why  other  subcortical  nuclei 
with  electrically  active  laminar  structures  cannot  be 
treated  in  the  same  manner.  Here  the  only  real  limita¬ 
tions  may  be  the  quality  of  signal-to-noise,  as  this  will 
be  degraded  for  very  deep  source  surfaces. 

It  is  obvious  that  the  next  step  will  be  to  implement 
the  extended  MNLS  procedures  using  actual  subjects 
and  images  of  their  brains.  Also,  clinical  trials  should 
be  attempted  where  there  are  confirmed  abnormal 
metabolic  levels  for  comparison  with  PET  ptxxedures, 
and  where  there  is  abnormal  field  power  in  some 
bandwidtbs  associated  with  epileptoid  phenomena. 
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